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Precise Half-Life Measurement of the 
Superallowed (3 + Emitter 26 Si 

V.E. IacobQ J.C. Hardy, A. BanuQ L. Chen, V. V. Golovko@ J. 
Goodwin, V. Horvat, N. NicaQ H.I. Park, L. Trache, and R.E. Tribble 

Cyclotron Institute at Texas A&M University 
(Dated: August 18, 2010) 

We have measured the half-life of the superallowed + — > + P + emitter 26 Si to be 2245.3(7) ms. We 
used pure sources of 26 Si and employed a high-efficiency gas counter, which was sensitive to positrons 
from both this nuclide and its daughter 26 Al m . The data were analyzed as a linked parent-daughter 
decay. To contribute meaningfully to any test of the unitarity of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix, the ft value of a superallowed transition must be determined to a precision of 0.1% 
or better. With a precision of 0.03% the present result is more than sufficient to be compatable with 
that requirement. Only the branching ratio now remains to be measured precisely before a ±0.1% 
ft value can be obtained for the superallowed transition from 26 Si. 

PACS numbers: 21.10.Tg,23.40.-s,27.30.+t 



I. INTRODUCTION 

The unitarity of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix is a fundamental requirement of the three- 
generation Standard Model. Currently, the most de- 
manding test available of CKM unitarity is the sum of 
squares of the experimentally determined elements of the 
matrix's top row [l|. The dominant term in this test is 
the up-down quark- mixing element, V U d, the most precise 
value of which is obtained through nuclear measurements 
of superallowed + — >Q + beta decays. To date, the mea- 
sured ft values for transitions from ten different nuclei 
are known to ~0.1% precision and three more to <0.3% 
So far, the superallowed transition from 26 Si has not 
been in either category, its /i-value precision being 0.8%, 
too large for it to contribute to the unitarity test. 

Since a superallowed + — transition involves only 
the vector current, its ft value relates to the vector cou- 
pling constant, G v - and, through it, to V u d - via the 
relationship [l[ 



Ft = ft{l + 5' R ){l + 5 NS -5 c ) = 



K 



2G 2 (1 



(1) 



where Ft is defined to be the "corrected" ft value and 
K/{hcf = 2tt 3 ft hi 2 /(m e c 2 ) 5 = 8120.2787(11) x 10~ 10 
GeV~ 4 s. There are four small correction terms: 5c 
is the isospin-symmetry-breaking correction; is the 
transition-independent part of the radiative correction; 
and the terms S' R and 5ns comprise the transition- 
dependent part of the radiative correction, the former 
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being a function only of the maximum positron energy 
and the atomic number, Z, of the daughter nucleus, while 
the latter, like 5c, depends in its evaluation on the details 
of nuclear structure. Both 5c and 5ns have been calcu- 
lated [2| with the best available shell-model wave func- 
tions, which are based on a wide range of spectroscopic 
data. They include those core orbitals that were deter- 
mined to be important based on measured spectroscopic 
factors in single-nucleon pick-up reactions; and they were 
further tuned to agree with measured binding energies, 
charge radii and coefficients of the isobaric multiplet mass 
equation. 

Although, in principle, the precise ft value for a 
single + — > + superallowed transition should be suffi- 
cient to determine V u d, that would leave the validity of 
these structure-dependent corrections without indepen- 
dent verification, and the derived value for G v could be 
in doubt. What has given credibility to the nuclear re- 
sult for G v is the fact that many superallowed transitions 
have been measured precisely and all give statistically 
identical results for Ft - and hence for G v . Since the 
uncorrected ft values actually scatter over a relatively 
wide range, it is the structure-dependent corrections that 
are responsible for bringing the Ft values into agreement 
with one another. 

Obviously, this is already a powerful experimental val- 
idation of the calculated corrections themselves, but it 
can be improved even further by precise measurements of 
additional transitions, especially ones calculated to have 
large correction terms. If the ft values measured for cases 
with large calculated corrections also turn into corrected 
Ft values that are consistent with the others, then this re- 
inforces the calculations' reliability for cases in which the 
corrections are smaller. The calculated correction terms 
for the superallowed transition from 26 Si give the result 
5c — 5ns = 0.65(3)% [2|. Although this value is not par- 
ticularly large compared to most of the well-measured 
cases, it is more than double the size of the correction 
for the superallowed transition from 26 Al m , which is its 
mirror transition. Experimental verification of this pre- 
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II. EXPERIMENT 



Overview 
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FIG. 1: Combined decay schemes of 26 Si and 26 Al m show- 
ing only those features of relevance to their superallowed /3 
decays. All energies are in keV and the Qec values shown 
are for the superallowed branch. The data are taken from 
Refs. 0,1. 



dieted mirror asymmetry would be a valuable test of the 
calculations. 

The ft value that characterizes any /3-transition de- 
pends on three measured quantities: the total transition 
energy, Qec, the half-life, ti/2, of the parent state and 
the branching ratio, R, for the particular transition of in- 
terest. The Q^c-value is required to determine the sta- 
tistical rate function, /, while the half-life and branching 
ratio combine to yield the partial half-life, t. The Qec 
value for 26 Si is already known sufficiently well Q to yield 
a value for / with 0.015% precision but, before the mea- 
surement reported here, the 26 Si half life was only known 
to 0.12% and its branching ratio to 0.8% 1]. Our new 
half-life for 26 Si has 0.03% precision and furthermore dis- 
agrees significantly with the measurement [ij that previ- 
ously dominated the world average. The result reported 
here represents our first step in bringing the precision of 
the ft value for this transition into the desired range of 
0.1%. The second step will be an improved branching 
ratio, a measurement that we will soon undertake. 

Like the decay of 34 Ar, whose half-life we have reported 
previously 26 Si /3 + decays to a daughter which is itself 
radioactive and is, in fact, another superallowed emitter. 
The combined decay schemes appear in Fig.[T] where it 
can be seen that the half-life of 26 Al m , the daughter, 
differs by only a factor of three from that of 26 Si. This 
adds complications to the experiment, which requires us 
to use the techniques that we developed previously for 
our study of 34 Ar. These will be described briefly in the 
following sections, but for further details the reader is 
referred to Ref. [5] . 



Precise half-life measurements need high-purity ra- 
dioactive beams, a requirement that is even more im- 
portant for cases where the daughter nucleus is also ra- 
dioactive with a similar half-life. We achieved this goal 
by using a production reaction with inverse kinemat- 
ics, p( 27 Al, 2n) 26 Si, and selecting the desired reaction 
product with the Momentum Achromat Recoil Separator 
(MARS) 0. Our primary beam of 30A-MeV 27 Al, which 
was produced by the superconducting cyclotron at Texas 
A&M University, impinged on a 1.6-atm hydrogen gas 
target cooled to liquid nitrogen temperature. The fully 
stripped ejectiles were then analyzed by MARS. Initially, 
working with a low-current primary beam, we inserted at 
the focal plane of MARS a 5 x 5 cm silicon telescope con- 
sisting of a 16-strip position-sensitive detector (PSD) 300 
fim. thick, backed by a 1-mm-thick detector. The tele- 
scope was used first for the identification of secondary 
reaction products, then for the control of the selection 
and focus of the desired species in the center of the beam 
line. This also gave us a clear indication of nearby reac- 
tion products that could potentially contribute as impu- 
rities to our selected beam. 

After the tuning and selection procedure, the PSD was 
dropped out of the way and the intensity of the primary 
beam was increased. With extraction slits at the MARS 
focal plane used to select 26 Si, the resulting 25.2-A-MeV 
radioactive beam was extracted into air through a 51-/zm- 
thick kapton window. This beam, typically of 2 x 10 4 ions 
per second, passed through a 0.3-mm thin BC-404 plas- 
tic scintillator, where the ions were counted, and then 
through a set of aluminum degraders, eventually being 
implanted in the 76-/im-thick aluminized mylar tape of a 
fast tape-transport system. The combination of m/q se- 
lectivity in MARS and range selectivity in the degraders 
led to implanted samples that were at least 99.8% pure. 

After the radioactive sample had been collected for a 
time interval of the order of the 26 Si half-life, the beam 
was turned off and the tape-transport system moved the 
sample in ~ 175 ms to a well-shielded location 90 cm away, 
stopping it in the center of a 4n proportional gas counter. 
The decay positrons were then detected for twenty half- 
lives (45 s), with signals from the gas counter being 
multiscaled into a 500-channel time spectrum. These 
collect-move-detect cycles were computer controlled and 
their timing was continuously monitored on-line. They 
were repeated, with a separate decay spectrum recorded 
for each, until the desired overall statistics had been 
achieved. In its shielded location, the gas counter had 
a background rate of about 0.5 counts/s, which was 3-4 
orders of magnitude lower than the initial count rate for 
each collected sample. For this experiment we accumu- 
lated data from 5,000 cycles split into 55 separate runs, 
which yielded a total of 2 x 10 8 counts. 
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B. Gas counter and electronics 

The gas counter we used is similar to ones we have 
used in previous precise half-life measurements 0, HHUl . 
It consists of two separate gas cells that, when assem- 
bled, have a 0.25-mm slot between them, through which 
the mylar transport-tape passes. Both cells were ma- 
chined from copper and each is equipped with anodes of 
13-/im-diameter gold-plated tungsten wire. Methane at 
just over one-atmosphere pressure is continuously flushed 
through both cells. Methane offers adequate gas gain for 
detecting positrons and is quite insensitive to 7 radia- 
tion. A Havar foil window, 3.7 cm in diameter and 1.5 
fim thick, hermetically seals each gas cell on the side fac- 
ing the tape. 

The electronic chain we used in the measurement was 
the same as that described in detail in Ref. The pre- 
amplified signal from the gas counter is first passed to a 
fast- filter amplifier with high gain (x500). At this high 
gain many of the pulses would saturate the amplifier so, 
to ensure that the amplifier recovers quickly, large pulses 
are clipped with a Schottky diode inserted after its first 
stage of amplification. The amplified and clipped pulses 
are then passed to a discriminator with very low thresh- 
old (150-250 mV). 

A 90 Sr/ 90 Y ft source is used to set up the detector 
system. This source has been specially prepared on a 
sample length of transport tape and is inserted into the 
gas counter in exactly the position that an on-line sam- 
ple, such as 26 Si, is situated. Our procedure is to record, 
at a fixed threshold setting, the counting rate from the 
discriminator as a function of the applied counter bias 
voltage. Initially, as the applied voltage is raised the 
count-rate also rises since the increasing gas gain leads 
to more primary ionizing ft events triggering the discrim- 
inator. However, at approximately 2600 volts - the exact 
value depends on the threshold setting - a "plateau" is 
reached, and the count rate remains nearly unchanged 
for the next 200-300 volts increase in the bias voltage. 
At higher voltages still, there is a second rapid rise in 
the count rate as spurious pulses increasingly trigger the 
discriminator. This behavior is well understood [12j and 
clearly demonstrates that, when operated in the plateau 
region, such detectors have essentially 100% efficiency. 
During our 26 Si measurement the detector was always 
operated in the plateau region. 

Since dead-time is a serious concern for half-life mea- 
surements, the discriminator output signals were split 
and sent to two fixed-width, non-retriggering and non- 
extending gate generators, which established different 
dominant dead times in the two separate streams, both 
of which were multiscaled. The time base for the mul- 
tiscalers was defined by a function generator, which is 
accurate to 0.01 ppm. Both gates also were continuously 
monitored during every run, thus giving us an on-line 
measure of the dead-time (±5 ns) during data collection. 
Note that even though the two gate generators were fed 
by the same data, the dead-time distortions of the un- 



derlying Poisson-distributed data are independent in the 
two cases; thus the two data streams allowed us to test 
that our dead-time corrected result was independent of 
the actual dead time of the circuit. 



C. Special precautions 

As the experiment was aimed at better than 0.1% pre- 
cision, many tests for systematic effects were made and 
special precautions taken during the measurements them- 
selves: 

• Every experiment was subdivided into many sep- 
arate runs, differing only in their particular com- 
bination of detection parameters: dominant dead- 
time, detector bias and discrimination threshold. 
We used combinations of four different dead times 
(3, 4, 6 and 8 fis), three discriminator thresholds 
(150, 200 and 250 mV) and four detector biases 
(2550, 2650, 2750 and 2850 V). The separate anal- 
ysis of each individual run allowed us to test for 
systematic effects that could contribute to the un- 
certainty in the final result. 

• Since each 26 Si decay produces an 26 Al m daughter 
that also decays, the relative activity of the two 
nuclides at the beginning of the detection period 
depends on the length of the collection period (and 
the tape-move time) that preceded it. We used four 
different collection times (1.0, 2.5, 3.0 and 4.5 s) to 
test for consistency. 

• The ratio of the parent to daughter activities also 
depends on the time-dependence of the rate at 
which 26 Si was accumulated in the tape during the 
collection period. The number of ions registered 
in the scintillator located just in front of the alu- 
minum degraders was recorded as a function of time 
with each cycle, and the results were used in our 
analysis. 

• The few weak impurities in the 26 Si beam have dif- 
ferent ranges in our degraders. Thus, any contribu- 
tion they might make to the total activity collected 
in the tape would depend on the depth at which 
the 26 Si beam is stopped in the tape. By using 
two different thicknesses of aluminum degrader we 
placed the 26 Si midway through the 76-/im tape for 
the first 30 runs and then near the rear surface for 
the remaining 25. Again we tested for consistent 
results. 

• The tape-transport system is quite consistent in 
placing the collected source within ±3 mm of the 
center of the detector, but it is a mechanical de- 
vice, and occasionally larger deviations occur. We 
separately recorded the number of implanted nu- 
clei detected in the scintillator during the collection 
period of each cycle, and the number of positrons 
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recorded in the gas counter during the subsequent 
count period. The ratio of the latter to the former 
is a sensitive measure of whether the source was 
seriously misplaced in the proportional counter. 

• We checked the composition of the beam exiting 
MARS once per day by re-inserting the PSD and 
ensuring that no changes had occurred. 

• Several background measurements were made, in 
which all conditions were identical to those of a 
normal run except that the tape motion was dis- 
abled. 

• In one run, we repeatedly collected activity for 16.5 
s and counted for 165 s in order to search for long- 
lived impurities. None was found. 

• It is important that the gas counter be operated in 
the "plateau" region: i.e. within the range of de- 
tector bias voltages in which the counting rate is 
nearly independent of voltage at a given discrim- 
inator setting. We measured this voltage plateau 
with a long-lived 90 Sr/ 90 Y source before and after 
our measurement, and once during it. In all cases 
we found the plateau slope to be <0.5% per 100 
V and did not observe any changes in the voltage 
boundaries of the plateau region. 

III. ANALYSIS AND RESULTS 

Before analyzing the data, we first removed any cycles 
that had fewer than 500 implanted 26 Si ions detected 
by the scintillator. These were the result of low - or 
no - primary beam current from the cyclotron during 
part or all of the collection period. Then we eliminated 
those cycles with an anomalously low ratio of recorded 
positrons to implanted silicon ions, which is indicative 
of faulty tape motion leading to a misplaced sample in 
the gas detector. Approximately 8% of the cycles were 
rejected for this reason. 

A. Parent-daughter connection 

For each run, we processed the accepted data by sum- 
ming the dead-time-corrected spectra from all its in- 
cluded cycles. We corrected for dead time using the pro- 
cedure outlined in Ref. @. The total time-decay spec- 
trum obtained from the combined runs is presented in 
Fig. [2] where we also show the separate contributions 
from the 26 Si parent and 26 Al m daughter. This break- 
down into components is based upon our final analysis 
and is presented here simply to illustrate how the parent- 
daughter decay curve, which combines two rather similar 
half-lives, tends to mask the parent half-life even though 
the parent activity dominates at the start of the counting 
period. 




Time [s] 

FIG. 2: Measured time-decay spectrum (solid line) for the 
total of all data obtained from the /3 + decay of 26 Si and its 
daughter 26 Al m . The dotted/dashed lines represent the de- 
rived 26 Si/ 26 Al m contributions. 

We can easily understand this situation by examin- 
ing the coupled decay equations for combined parent- 
daughter decays. The combined 26 Si and 26 Al m activity 
yields a total rate for detected positrons of 

Am = C ie - Alt + C 2 e- A2t , (2) 

with 

Ci = iWr ( * T \-) 
\e 2 Ai-A 2 / 

<* = JW»(f + Y^)' (3) 

where t is the time elapsed after the end of the collect 
period; Nip are the numbers of 26 Si and 26 Al m nuclei 
present in the sample at t = 0; ei^ are the experimental 
efficiencies for detecting the positrons from the respective 
decays; and Ai^ are the corresponding decay constants. 
Note that when Ai = 2A2 (and t\ — £2) the coefficient 
Ci vanishes, leaving a single exponential term having the 
decay constant of the daughter. The half-lives of 26 Si 
and 26 Al m are actually related by a factor of 2.8, close 
enough to 2 that, for our measurements, the coefficient 
C\ was slightly smaller than C2, leaving the daughter 
to dominate the decay curve (see Fig. [5]). This imposes 
a real limitation on the precision that can be obtained 
from a conventional fit to the data: even with A2 fixed 
at its known value, C±, C2 and Ai (as well as the con- 
stant background) must all be determined independently, 
which leads to Ai having a rather large uncertainty. 

We can overcome this limitation by fixing the ratio 
C2/C1 so as to reduce the number of adjustable parame- 
ters in the fit to three (including background), as we did 
for our measurement of the 34 Ar half-life [5j. In practice 
this means that we must establish two ratios, N2/N1 and 
ei/e2, from the experimental parameters. In this section 
we deal with the former; in the next section, with the 
latter. 
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FIG. 3: Typical time-profile of the collected 26 Si beam. The 
initial drop in intensity is generated by the change in gas 
density as the primary beam heats the gas around its path. 
A fan located inside the gas-target ensures a rapid transition 
to stable conditions. 

No 26 Al ions were present in the 26 Si sample collected 
in our tape (see Sec. lIII Cj) . Thus, if the sample col- 
lection rate were exactly constant, we could determine 
N2/N1 from a simple calculation of the production of 
26 Al m (via 26 Si decay) over the precisely known time of 
the collection period. However, we measured the actual 
number of 26 Si ions as a function of time with the scin- 
tillator at the exit of MARS and determined that there 
was a slightly higher rate at the very beginning of the 
collection period (see Fig. [3]). Our cryo-target was a gas, 
and the primary beam heats that gas around its path 
through the target, thus initially generating a local drop 
in gas density; the transition to steady conditions was 
hastened by a fan that continuously circulated the gas 
in the target cell. We also found that the size of the 
variation in the radioactive-beam intensity at the begin- 
ning of a cycle depended on the primary beam intensity, 
potentially changing the beam-profile from one cycle to 
another. Nevertheless, with the collection time-profile 
measured and recorded for each cycle, we could perform 
a numerical integration over the measured 26 Si accumu- 
lation to calculate the decay-production of 26 Al m and the 
corresponding N2/N1 ratio accurately for each run. 

B. Parent-daughter relative efficiencies 

The ratio e\/t2 can also be established independently. 
In our experiment, we detected positrons from the de- 
cays of both 26 Si and 26 Al m with a very low threshold 
and nearly 100% overall efficiency, so superficially one 
might conclude that ei/e2 = 1 and that the efficiency 
ratio can consequently be ignored. However, the end- 
point energy of the (3 + spectrum corresponding to the 
26 Si superallowed transition is 3819 keV, while that for 
the equivalent transition from 26 Al m is 3211 keV. The 
608-keV difference between them means that the shapes 
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FIG. 4: (a) Calculated /3 + spectra for 26 Si (solid curve) and 
26 Al m (dashed curve). The former includes the Gamow- Teller 
branches to 1 + states as well as the superallowed branch to 
the + isomer; the latter is pure superallowed decay (see Fig.fj] 
for both decay schemes), (b) System efficiency for detecting 
positrons due to the effects of the aluminized Mylar tape and 
the Havar windows of the detector gas cells. The curve is the 
result of a Monte Carlo calculation. 



of the two spectra are slightly different from one another. 
Furthermore, the Gamow- Teller branches that occur only 
for the 26 Si decay contribute to the differences as well. 
This is illustrated by the calculated (3 + spectra shown in 
Fig.Ha). 

Because we operate our gas detector in the plateau 
region (see Sec. Ill Bp . we can be confident that our elec- 
tronic detection threshold is too low for losses from that 
source to have any effect on our overall detection efficien- 
cies. However, the aluminized Mylar tape (half-thickness, 
38 /xm) and the Havar window of each detector gas cell 
(1.5 /im thick) stop the lowest-energy positrons, thus pre- 
venting them from reaching the active volume of the de- 
tector. In effect, this imposes a low-energy cut-off and, 
since the parent and daughter j3 spectra have slightly 
different shapes at low energies, the fraction of positrons 
lost in one case is slightly different from that in the other. 
For a pair of decays like ours, where the average positron 
energy is greater for the parent than it is for the daugh- 
ter, the ratio is always t\/e2 > 1. 

Since the decay positrons are emitted isotropically, 
their paths through the tape and window cover a range 
of lengths, resulting in a low-energy cut-off that is not 
sharp. Nevertheless, for any given cut-off energy the 
effect on the efficiency ratio can readily be calculated 
from the known /3-spectrum shapes. Using the Monte 
Carlo code EGSnrc (version V4-r2-3-0) [l3| , in which we 
modeled our exact tape/ window/detector geometry, we 
obtained the overall system efficiency as a function of 
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positron energy, the result being shown in Fig.[4|b). Note 
that the code was only required to calculate the losses 
due to the tape and window. We have extensively tested 
the accuracy of three Monte Carlo codes - EGSnrc, 
Geant4 and PENELOPE - in fitting experimental results 
from a thin scintillation detector for low energy conver- 
sion electrons as well as /3-decay spectra. We found that 
EGSnrc offers the best combination: it agrees well with 
experiment and it operates most efficiently. 

The Monte Carlo result was then integrated with each 
of the two spectra in Fig.[4|a) to obtain our overall effi- 
ciencies for detecting the parent and the daughter ac- 
tivities. From these, we derived the efficiency ratio, 
ei/e2 = 1.00143(25), where we have assigned an uncer- 
tainty that encompasses a ±10% relative uncertainty in 
the calculated ranges. This uncertainty was based on an 
assessment of how much the ratio would change if the 
source were not located at exactly the center of the de- 
tector, thus slightly changing the losses in the tape and 
window. However, by eliminating all cycles in which the 
ratio of recorded positrons to implanted silicon ions was 
anomalously low, we ensured that all analyzed cycles cor- 
responded to central or nearly central source locations. 
The uncertainty we have assigned to the efficiency ratio 
is a very generous allowance for the range of locations 
allowed in the accepted cycles. The value we obtained 
for ei/e2 and its uncertainty were subsequently used in 
the analysis of all runs. 



C. Sample impurities 

From the position-sensitive detectors that we inserted 
into the MARS focal plane before the measurement and 
periodically during it, we could identify and monitor any 
reaction products that might potentially contribute as 
impurities to our selected 26 Si beam. Only three were 
detectable and all were extremely weak: 25 Al (0.02% of 
the 26 Si intensity), 24 Al (0.14%) and 23 Mg (0.02%). How- 
ever, even these small percentages must be reduced when 
we consider what was actually retained in the aluminized 
Mylar tape that transported the collected activity to our 
detector. As the impurity ions passed through our alu- 
minum degraders, each impurity experienced a different 
energy loss from the others and from 26 Si. The result is 
illustrated in Fig. [3] for the two thicknesses of degraders 
employed during the measurement (see Sec.Ql]). 

Figure [S] presents the results of calculations based on 
the SRIM code 115]. Before the main measurements be- 
gan, we recorded the collected 26 Si activity as a function 
of aluminum degrader thickness and were thus able to 
determine experimentally the precise thickness required 
to center the 26 Si deposit in the tape. The value ob- 
tained was very close to that predicted by SRIM, which 
gives confidence that the calculations can be relied on to 
determine the spatial distributions of the impurities rel- 
ative to that of 26 Si. From Fig. [5] it is evident that only 
10% of the 24 Al and 50% of the 23 Mg was collected in 
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FIG. 5: Illustrations of calculated implantation profiles in 
Mylar for the 26 Si beam and those impurities with similar 
ranges. All beams enter from the left. The top illustration 
(a) gives the profile after the beams have passed through 159 
lira of aluminum, while the bottom one (b) gives the result 
after 171 fira of aluminum. The shaded region in both illus- 
trations corresponds to the actual thickness of our collection 
tape. Those ions within the shaded region are collected in our 
sample; all others are not. 

the tape when the 26 Si was centered in the tape (top il- 
lustration) and none at all when the 26 Si was placed near 
the back of the tape (bottom). Consequently, relative to 
the collected 26 Si activity, the collected 24 Al and 23 Mg 
activities were both approximately 0.01% when the 26 Si 
was centered; they were zero when it was near the back. 
The only remaining impurity, 25 Al, was fully collected at 
both degrader settings, but its relative intensity was only 
0.02% to start with. The effect of these three impurities 
was only barely significant, but we included them with 
their appropriate relative intensities when we fitted the 
collected decay spectra. We also incorporated a ±30% 
relative uncertainty on the intensities used. 

D. Time decay analysis 

We fitted the data from each of the 55 runs sepa- 
rately, incorporating four components: 26 Si; its daugh- 
ter, 26 Al m ; the weak impurities, 25 Al, 24 Al and 23 Mg; 
and a constant background. The half-life of 26 Al™ 1 was 
fixed at its known value of 6345.0(19) ms [l|. Its ini- 
tial activity relative to that of 26 Si was set for each run 
to the value obtained from numerical integration of the 
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FIG. 6: Test for possible systematic bias in the 26 Si half- life 
measurement due to discriminator threshold or detector volt- 
age. Open/grey /black symbols represent the three discrimi- 
nator settings, 150 mV/200 mV/250 mV; the four detector bi- 
ases, 2550 V, 2650 V, 2750 V and 2850 V are represented by 
the symbol shapes A, □, O and V, respectively. The average 
value for the half-life is 2245.26(51) ms (statistical uncertainty 
only) with x 2 /ndf = 70.6/54. The average value appears as 
the solid line, with dashed lines as uncertainty limits. 
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FIG. 7: Test for possible systematic bias in the 26 Si half-life 
measurement caused by undetected short-lived impurities or 
by rate-dependent counting losses. Each point is the result of 
a separate fit to the data; the abscissa for each point repre- 
sents the time period at the beginning of the counting cycle 
for which the data were omitted from that particular fit. The 
solid and dashed lines correspond to the average half-life value 
and uncertainties given in Fig. [6] 



measured time-profile for the collected 26 Si beam in that 
run (see Sec. MI Al and Fig. [3]); and the efficiency ratio, 
ei/e 2 , was fixed at the value established in Sec. lIIIBl 
The half-lives of 25 Al, 24 Al and 23 Mg were taken to be 
7.182(12) s HI, 2.053(4) s @ and 11.324(10) s [H respec- 
tively. Their initial relative activities were obtained from 
the measurement and calculations described in Sec. lIII Cl 

Since each run was obtained with a different combina- 
tion of detection settings, we could use the individually 
fitted 26 Si half-lives to test for any systematic depen- 
dence on those settings. As seen in Fig. [51 the half-life 
results show no systematic dependence on detector bias 
or discriminator threshold. Although not illustrated, the 
results were also found to be independent of both the im- 
posed circuit dead time and the length of time for which 
the sample was collected. As a final systematic check, 
in this case for the possible presence of short-lived impu- 
rities or other possible count-rate dependent effects, we 
removed data from the first 0.9 s of the counting period 
in each measurement and refitted the remainder; then we 
repeated the procedure, removing the first 1.8 s, 2.7 s and 
so on. From Fig. [7] it can be seen that, within statistics, 
the derived half-life was also stable against these changes. 

With these possible systematic effects eliminated as 
significant factors, we can combine all 55 runs to obtain a 
value for the 26 Si half-life of t 1/2 ( 26 Si) = 2245.26ms, with 
a statistical uncertainty of ±0.51 ms. The normalized \ 2 
of this average is 1.31 and the statistical uncertainty has 
been multiplied by the square root of this number. The 
result itself represents a self-consistent analysis of about 
200 million combined 26 Si and 26 Al m decay events. 

This analysis has depended upon the source being cen- 
trally located in the gas counter for each cycle since 
serious misplacement would have resulted in a differ- 



ent parent-daughter relative efficiency for that cycle (see 
Sec. lIII Bj) . We ensured centrality by rejecting cycles 
that had anomalously low ratios of recorded positrons to 
implanted silicon ions. For each run we included only 
those cycles with ratios between 91 and 100% of the 
maximum value obtained for that run. To test whether 
our result is in any way sensitive to this choice, we re- 
peated the full analysis for subgroups of the cycles cor- 
responding to ratios between 98-100%, 96-98%, 94-96% 
and so on. Although having larger uncertainties, the half- 
lives obtained from those subgroups within our selected 
band of 91-100% were all statistically consistent with our 
quoted result of 2245.3(5) ms. More importantly, the 
nearest group of rejected cycles - those with ratios in the 
range 88-90% - were also consistent, yielding a value of 
2246.7(30) ms. Clearly, our selection criterion for accept- 
ing cycles was a conservative one that does not adversely 
affect the result. 

It is also interesting to compare this result with the 
half-life value derived from a fit to the data that does 
not impose a fixed link between the parent and daugh- 
ter activities. We described in Secs. lIII Al and MI Bl how 
we could reduce by one the number of free parameters 
in the fit by independently determining the ratios N2/N1 
and ei/e2- This improved the statistical uncertainty in 
the fitted result but did introduce an additional uncer- 
tainty associated with the ratio ei/e2 (see the second line 
in Tablefl]) . If we ignore the parent-daughter link and fit 
the data with four variable parameters (including back- 
ground), we obtain a half-life result of 2243.2(22) ms, 
which is statistically consistent with the value presented 
in Table|T]but has a thrice larger uncertainty. 



8 



E. Uncertainty budget 

There are other contributions, of course, to our fi- 
nal uncertainty beyond the contribution from counting 
statistics. We itemize them all in Table Q] Counting 
statistics is the largest contributor to the overall uncer- 
tainty, but the contribution associated with the different 
efficiencies for detecting parent and daughter activities 
is a close second. Less important contributors were the 
uncertainties in the 26 Al m half-life, the dominant circuit 
dead-time and the weak sample impurities. Our final 
result for the 26 Si half-life is 2245.3(7) ms, in which sta- 
tistical and systematic uncertainties have been combined. 



F. Comparison with previous results 

Three previous measurements of the 26 Si half-life are 
listed in the most recent review [l[: 2210(2f) ms [l7j j . 
2240(10) ms [H and 2228.3(27) ms Q. The first two 
of these results were recorded more than 30 years ago 
and have uncertainties larger by more than a factor of 
10 than our present result; one of them agrees with our 
result, while the other lies one-and-a-half of its standard 
deviations away. This can hardly be viewed as a matter 
of concern. However, the last of the previous measure- 
ments was published only two years ago, has a quoted 
uncertainty that is only four times ours and differs from 
our result by more than six times that uncertainty. This 
discrepancy has to be addressed. 

In their measurement, Matea et al. [4] employed a pu- 
rified 26 Si beam from the JYFLTRAP Penning-trap fa- 
cility and implanted it in a movable 100-/^m-thick Mylar 
tape located at the center of a hollow cylinder whose 
2-mm-thick walls were made from plastic scintillator op- 
tically coupled to two photomultipliers. After sample im- 
plantation, the decay positrons were counted in the plas- 
tic scintillator, after which the tape removed the sample 
and the cycle was repeated. The experimental arrange- 
ment is pictured in Ref. [l9j]; there it can be seen that 
the tape enters and exits the scintillator cylinder via the 
same opening and passes over a roller, which is located 
within the cylinder. 

The method Matea et al. used to analyze their data 



TABLE I: Error budget for the 26 Si half- life measurement. 



Source 



Uncertainty (ms) 



statistics 

efficiency ratio, ei / £2 
t 1/2 ( 26 Al m ) 
dead time 
sample impurities 

Total 

26 Si half-life result (ms) 



0.51 
0.37 
0.16 
0.07 
0.04 

0.66 

2245.3(7) 



was similar to ours in that they fixed the activity of 
26 AP™ based on its calculated production from the de- 
cay of an initially pure 26 Si sample. Unfortunately they 
did not take account of any detection-efficiency differ- 
ence between the parent and daughter activities [20| (see 
Sec. lIII Bj) . Furthermore, their experimental conditions 
make this a much more serious problem for them than 
it was for us. First, both their thicker tape and the 
roller they used to allow that tape to move introduced 
a more significant and complicated low-energy cut-off 
due to positrons ranging out in those materials. Second, 
with a plastic scintillator and photomultiplier tubes gen- 
erating their signals, they certainly needed to set a non- 
negligible low-energy electronic threshold. Third, their 
scintillator was, to some extent, sensitive to 7 rays, which 
would slightly favor detection of 26 Si, which produces /?- 
delayed 7 rays, while 26 Al'™ does not. 

Although it is obviously impossible for us to model in 
detail someone else's experimental apparatus, we approx- 
imated their arrangement in the EGSnrc Monte Carlo 
code [IH, the same code we used to help determine the 
efficiency ratio for our measurement (see Sec. lIII Bj) . We 
found that the results were very sensitive to just those 
experimental details that we could not define precisely. 
Nevertheless, we determined that the value of ei/e2 could 
easily reach 1.005 or possibly larger. If this ratio had 
been applied to their data, it would have increased their 
measured half-life by ~0.4%, bringing it to 2237 ms, 
approximately half way to our result, and undoubtedly 
with a much larger uncertainty than the 2.7 ms origi- 
nally claimed. Clearly this is only an indicative not a 
definitive calculation, since the exact correction and its 
application to the data requires information only avail- 
able to the original authors. For now, though, we believe 
that it reveals a serious omission in the original analysis 
by Matea et al. [|| , and that their published result should 
simply be discounted. 



IV. CONCLUSIONS 

We report here the first measurement of the half-life 
of the superallowed j3 + emitter 26 Si obtained with a pre- 
cision of better than 0.1%. Since 26 Si and its daughter 
26 Al m have comparable half-lives, this measurement re- 
quired us to use the technique we developed for the mea- 
surement of the 34 Ar half-life in which we link the 
parent and daughter decays based on a precise knowledge 
of the rate of deposition of the 26 Si source. We also had 
to make small but important corrections for detection- 
efficiency differences between the parent and daughter 
/9 + decays. Possible sources of error were carefully inves- 
tigated and an error budget established. 

Our new precise result for the 26 Si half-life, 2245.3(7) 
ms, is considerably different from the average value 
quoted in the most recent survey of world data for su- 
perallowed + ^ + /3-decay transitions [l[. However, the 
average there is dominated by a measurement Q that 
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we now argue is flawed and should be discarded, since it 
leaves out the correction for parent-daughter detection- 
efficiency differences. 

With our half-life result for 26 Si decay having 0.03% 
precision, and the Q ec value for its superallowed branch 
being recently determined to 0.0025% 3], the Ft value 
for the branch can in principle be determined to sub- 
0.1% precision if the branching ratio can be measured 
to that level. This is a difficult but potentially achiev- 
able goal and, if accomplished, would bring 26 Si decay 
into the same category of precision as the best known 
superallowed transitions. 

This would allow a particularly interesting comparison 
with its mirror, 26 Al m decay. If calculated with Saxon- 
Woods radial wave functions (H, the nuclear-structure- 
dependent correction, {8c — 5ns), for the 26 Al m decay is 
0.305(27)% [J, the smallest value for any superallowed 
transition, while the equivalent correction term for the 



26 Si decay is 0.650(34)%. The difference between them, 
0.345(43)%, is reduced considerably, to 0.145(81)%, if 
Hartree-Fock wave functions are used instead [l| ■ An ex- 
perimental determination of the difference between these 
mirror transitions with sub-0.1% precision would discrim- 
inate between these theoretical approaches and could po- 
tentially reduce the theoretical uncertainties that affect 
the ultimate determination of V u d and the unitarity of 
the CKM matrix. 
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